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Abstract 

We calculate pionic contribution to the relativistic Fermi Liquid parameters (RFLPs) using Chiral Effective Lagrangian. The RFLPs 
^SJ so determined are then used to calculate chemical potential, exchange and nuclear symmetry energies due to t:N interaction. We 

also evaluate two loop ring diagrams involving a, uj and n meson exchanges and compare results with what one obtains from the 
' ^ 'relativistic Fermi Liquid theory (RFLT). 
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Fermi liquid theory (FLT) provides us with one of the 
most important theoretical scheme to study the properties 
of strongly interacting Fermi systems involving low lying 
excitations near the Fermi surface pQ . Although developed 
originally in the context of studying the properties of 3 He, 
it has widespread applications in other disciplines of many 
body physics like superconductivity, super fluidity, nuclear 
and neutron star matter etc |2|3j . 

In nuclear physics FLT was first extended and used by 
Migdal [4] to study the properties of both unbound nuclear 
matter and Finite nuclei [5] . FLT also provides theoretical 
foundation for the nuclear shell model [5] as well as nuclear 
dynamics of low energy excitations |2|6j . Particularly, ref. [7] 
reveals the connection between Landau, Brueckner-Bethe 
and Migdal theories, ref. [819] , on the other hand, calculates 
the Migdal parameters using one-boson-exchange models of 
the nuclear force and shows how these parameters are mod- 
ified if nuclear matter is considered within the context of 
Relativistic Brueckner-Hartree-Fock (RBHF) model. Some 
of the recent work that examines Fermi-liquid properties of 
hadronic matter also incorporates Brown-Rho (BR) scal- 
ing which is very important for the study of the properties 
of hadrons in dense nuclear matter [10) 1 lj . 

Most of the earlier nuclear matter calculations that in- 
volved Landau theory were done in a non-relativistic frame- 
work. The relativistic extension of the FLT was first de- 



Email addresses: kausik. paiasaha. ac . in (Kausik Pal), 
subhrajyoti .biswas<3saha. ac . in (Subhrajyoti Biswas), 
abhee .dm9saha. ac . in (Abhee K. Dutt-Mazumder). 



veloped by Baym and Chin [3] in the context of studying 
the properties of dense nuclear matter (DNM). In [3] the 
authors invoked Walecka model to calculate various inter- 
action parameters {f PP ' ) but did not consider mean fields 
(MF) for the a and to meson i. e. there the FLPs are calcu- 
lated perturbatively. 

Later Matsui revisited the problem in [12] where one 
starts from the expression of energy density in presence of 
scalar and vector meson MF and takes functional deriva- 
tives to determine the FLPs. The results are qualitatively 
different than the perturbative results as may be seen from 
[3 12 . A comparison of relativistic and non-relativistic cal- 
culations have been made in [8|9j which also discusses how 
the FLPs are modified in presence of the a and u> MF and 
contrast those with perturbative results. 

Besides a and u> meson, ref. [12] also includes the p and 
7T meson and the model adopted was originally proposed 
by Serot to incorporate pion into the Walecka model. It 
is to be noted, however, that the FLPs presented in [T^] 
are independent of tt meson. This is because 7r, being a 
pseudoscalar, fails to contribute at the MF level. Hence to 
estimate the pionic contribution to FLPs it is necessary to 
go beyond MF formalism. 

It is to be noted that, to our knowledge, such relativistic 
calculations including tt exchange does not exist, despite 
the fact that pion has a special status in nuclear physics as it 
is responsible for the spin-isospin dependent long range nu- 
clear force. Furthermore, there are various non-relativistic 
calculations including the celebrated work of Migdal which 
shows that the pionic contribution to FLPs are important 



Preprint submitted to Elsevier 



24 March 2008 



[10111113] and most dominant for low energy excitations. It 
might be mentioned here that [TT] discusses how can one 
incorporate relativistic corrections to F™ in the static po- 
tential model calculation. We, however, take the approach 
of [3] where all the fields are treated relativistically. 

The other major departure of the present work from 
ref. |3U2j resides in the choice of model for the description 
of the many body nuclear system. The straightforward in- 
corporation 7T meson into the Walecka model as was done 
in ref.pT] has serious difficulties. In particular, at the MF 
level it gives rise to tachyonic mode in matter at densities 
as low as O.lpo, where po is the nuclear saturation den- 
sity |15] . Inclusion of exchange diagrams removes such un- 
physical mode but makes the effective mass unrealistically 
large. Replacing pseudoscalar coupling with the pseudovec- 
tor interaction avoids these difficulties, however, it turns 
the theory non-renormalizable [IB]. The theoretical chal- 
lange therefore was then to construct a model with irN 
interaction which preserves renormalizibility of the theory 
and at the same time yields realistic results for the pion 
dispersion relations in matter. This was accomplished in 
[IB] , But this model was also not found to be trouble free, 
particularly it had several shortcomings in describing irN 
dynamics in matter which we do not discuss and refer the 
reader to |17I18I19I2()I2TT22] . 

These apart, Walecka model itself has several problems 
in relation to convergence which forbids systematic expan- 
sion scheme to perform any perturbative calculations. This 
was first exposed in rcf . [23 . Since then a lot of theoretical 
efforts have been directed to circumvent these problems. 
The most recent model that cure all these maladies and 
provides us with a systematic scheme to study the dense 
nuclear system is provided by Chiral Effective Field the- 
ory (xEFT) [21123124125] . The xEFT, apart from er and u 
mesons, also includes pion and therefore is best suited for 
the present purpose. 

In this letter we estimate contributions of pion exchange 
to the FLPs within the framrework of RFLT and subse- 
quently use the parameters so determined to calculate var- 
ious quantities like pionic contribution to the chemical po- 
tential, energy density, symmetry energy (/?) etc. For com- 
pleteness and direct comparison with the two loop results 
we also calculate here the exchange energy due to the in- 
teraction mediated by the er, lo and n mesons. The pionic 
contribution to the effective chemical potential, as we shall 
see, is significantly large. This, therefore, might alter the f3 
equilibrium condition than what one obtains in MF calcu- 
lations [12] . This we plan to study in a separate work [26] . 

Now we quickly outline the formalism. In FLT total en- 
ergy E of an interacting system is the functional of occupa- 
tion number n p of the quasi-particle states of momentum p. 
The excitation of the system is equivalent to the change of 
occupation number by an amount Sn p . The corresponding 
energy of the system is given by [213] , 
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where E is the ground state energy and s is the spin index, 
and the quasi-particle energy can be written as, 
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where e ps is the non-interacting single particle energy. It is 
to be remembered, that, although the interaction (f ps ,p>s') 
between the quasiparticles is not small, the problem is 
greatly simplified because it is sufficient to consider only 
pair collisions between the quasiparticles [5] . 

Since quasi-particles are well defined only near the Fermi 
surface, one assumes 
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Then LPs fis are defined by the Legendre expansion of 
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where 9 is the angle between p and p' , both taken to be on 
the Fermi surface, and the integration is over all directions 
of p. We restrict ourselves for I < 1 i.e. f and /i, as higher 
/ contribution decreases rapidly [3112] . 

Now the Landau Fermi liquid interaction f ps ,p's' is re- 
lated to the two particle forward scattering amplitude via 
0, 

Jps,p's' — n Atp^p's', 

£ p t p' 



(5) 



where the Lorentz invariant matrix M ps .p's' consists of the 
usual direct and exchange amplitude, which may be eval- 
uated directly from the relevant Feynman diagrams. The 
spin averaged scattering amplitude (f pp ') is given by [3], 
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The dimensionless LPs are Fi = N(0)fi, where N(0) is 
the density of states at the Fermi surface defined as [T2] , 
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Here g s , gi are the spin and isospin degeneracy factor re- 
spectively. 

By retaining only the lowest order terms in the pion 
fields, one obtains the following Lagrangian from the chi- 
rally invariant Lagrangian [24125] : 
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whereof = d^us u — d^ui^, tt_ = ^(r-^^) and f is the isospin 
index for pion. Here \& is the nucleon field and $o-, oj^ and 
$ w are the meson fields (isoscalar-scalar, isoscalar-vector 
and isovector-pseudoscalar respectively). The terms Cnl 
and SC contain the non-linear and counterms respectively 
(for explicit expressions see [25]). 

Now due to presence of pion fields in the chiral La- 
grangian we have component in the interaction which acts 
on the isospin fluctuation. One can derive the isospin de- 
pendent quasiparticle interaction along the line of ref.[3]. 
For pions, as mentioned before, the direct term vanishes 
and it is only the exchange diagram that contributes to the 
interaction parameter : 
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where g\ = f .5876 , f„ = 93MeV and = 139 MeV[ 
The factor 3 arises since 
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It is to be noted that, in the massless limit for a and u> 
meson, f^ x and /-f x diverge as shown in [5] , in contrast for 
pion, even in the massless limit, these are finite. This is due 
to the presence of (1 — cos#) in the numerator of Eq.© 
unlike a and u> meson. 

The dimensionless LPs can be determined by equating 
the equation Fq = A^(0)/^ and F{ = N(0)f[, where N(0) is 
the density of states at the Fermi surface defined in Eq. ([7]) . 
Thus the dimensionless parameters arc 
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(17) 



Ti ■ f 2 — 2(ti + t 2+ + Ti_r 2 -) + T10T20, 



(10) 



where r + and r_ are the isospin raising and lowering oper- 
ators, respectively, and tq represents the third component 

of f [13]. 

The effective nucleon mass M* is determined self- 
consistently from the following equation [12] , 
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Using Eq. ^ and Eq. © we can derive isovector LPs f' Q and 
fl, 
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Using Eqs.([T2]) and ([13]) we find that 
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It is this combination i.e. fo — 5/1, which appears in the 
calculation of chemical potential and other relevant quan- 
tities. For the massless pion, Eq. ([T4]) turns out to be finite, 




Fig. 1. Dimensionless Landau parameters in symmetric nuclear mat- 
ter for pion exchange in rclativistic theory. Solid and dashed line 
represent F^ x and F£ x respectively. 

In Fig.fl]) we show the density dependent of F$ x and Ff x 
due to pionic interaction. Numerically at nuclear matter 
density (p = 0.148ml- 1 ), = 1.35 and F' x ex ^ = 

—0.40. This might be contrasted with what one obtains in 
various other models. For example, at p = po [10] and [13] 
report F x x '* = —0.45 and —0.35 respectively. 

We now proceed to calculate the chemical potential due 
to the exchange terms denoted by p ex . As in ref.[3] we have 
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To calculate p, it is sufficient to let p = Sf in the right hand 
side of Eq. ([T8]) . With the constant of integration adjusted 
so that at high density pf — £/, Eq. (|18p upon integration 
together with Eq. ([T4]) yield 
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Table 1 

Dimensionless Landau parameters and chemical potential at p = po. 



Meson 


F 


Fi 


H ex (MeV) 


(T 


-5.04 


0.875 


731.89 


U> 


5.44 


-0.93 


501.82 


7T 


1.35 


-0.40 


635.78 
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where x = pf/M* and = m^/M*. 

The calculations of LPs and exchange chemical potential 
for other mesons is straightforward. However, for brevity, 
we do not present corresponding expressions for a and uj 
mesons but quote their numerical values in Table([l]). The 
numbers cited above are relevant for normal nuclear matter 
density po = 0.148fm -1 . For the coupling constants we 
adopt the same parameter set as designated by MOA in 

Interestingly, individual contribution to LPs of a and oj 
meson are large while sum of their contribution to F$ ot is 
small due to the sensitive cancelation of Ffi and Ftf as can 
be seen from Tablc([T]). Such a cancellation is responsible for 
the nuclear saturation dynamics |8l9j . Numerically, Fq +uj 
is three times smaller than Fq. 

Once the is determined, one can readily calculate its 
contribution to the energy density |3|27j . 
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where r\ = \fl + x 2 and 



/- = -2ylV^f- tan" 1 (^EK) 
For the massless pion this reads as 
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The contribution arising from pion exchange from the 
direct evaluation of Fig|2](c) reads as [25] . 
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Fig. 2. Two-loop contributions to the nuclear matter energy density. 
The solid line represents the baryon propagator, cr, u> and n mesons 
are denoted by dotted, wavy and dashed line respectively. 

Table 2 

Exchange energy in MeV from two loop diagram and FLT at p = po. 



Meson 




E ex 




Two-loop 


FLT 


a 


42.18 


40.48 


UJ 


-28.17 


-23.41 


7T 


13.66 


12.49 
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Similarly from FigO-((a),(b)) one can determine exchange 
energy due to a and u> meson interaction [25]. In Table© 
we compare the exchange energy results obtained from the 
direct evaluation of two loop diagrams with those calcu- 
lated from RFLPs. 

Knowing the isovector LPs, to which here only the pion 
contributes, one can calculate nuclear symmetry energy. 
The symmetry energy is defined as the difference of energy 
between the neutron matter and symmetric nuclear matter 
is given by the following expression |12l28j 



1 d 2 E, 
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In terms of LPs, the symmetry energy can be expressed as 
[12 28 



(25) 



Using Eq. (fT6| and Eq. (|25|) , numerically at saturation den- 
sity (p = po) we obtain (3 = 44.36 MeV. 

In this letter, we calculate RFLPs within the framework 
of RFLT. For the description of dense nuclear system %EFT 
is invoked. Although our main focus was to estimate the 
contribution of pions to RFLPs, for comparison and com- 
pleteness we also present results for the a and ui meson. It 
is seen that the pionic contribution to the FLPs are signifi- 
cantly larger compared to the combined contributions of a 
and lu meson. Thus any realistic relativistic calculation for 
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the FLPs should include ir meson which necessarily implies 
going beyond the MF calculations. The LPs what we de- 
termine here are subsequently used to calculate exchange 
and symmetry energy of the system. Finally we evaluate 
two loop ring diagrams with the same set of interaction pa- 
rameters and show that the numerical results are consistent 
with those obtained from the FLT. 
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